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Abstract 

We study thermodynamic properties of an Ising model of a ferromagnetic nanoscopic pyramid 
deposited onto a ferromagnetic bulk substrate. The influence of the interaction between the pyra- 
mid and the substrate is calculated in terms of the equilibrium reduced-state (density) operator 
used for description of thermodynamic properties of nanoscopic systems. The spatial distribution 
of the fluctuations of molecular field and magnetization in the nanoscopic pyramid is obtained 
in the Gaussian fluctuations approximation. Experimental consequences for the magnetic force 
measurements are briefly discussed. 



PACS numbers: 75.75.+a 



1. INTRODUCTION 



During the past two decades we have witnessed significant advances in the abihty to syn- 
thesize nanoscale structures as well as development of novel experimental method allowing 
exploration of their physical properties. |l| This is exciting for two reasons. Firstly, new 
forms of matter with no counterpart in nature and revealing unique physical properties have 
been fabricated. Secondly, we have now realized nanostructures that open new avenues for 
development of very small devices. 

Nanoscopic magnetic systems[2| have always been very attractive from the theoretical 
point of view. On the other hand, studies of nanoscopic systems are characterized by a close 
coupling between theory and experiment because of rapidly increasing number of experi- 
mental works on real nanoscopic materials. [l| These materials often correspond remarkably 
close to certain idealized spin models and are great challenges not only for physicists, but 
also for chemists and engineers. 

There is a great interest in producing two-dimensional arrays of magnetic nanodots 
which may serve e.g. as magnetic recording media. Due to small dimensions of these 
particles quantum-mechanical effects cannot be ignored. Effects observed in quantum dots 
are relevant also in quantum computing and spin electronics. It was shown recently that 
nanoscopic structures in the form of piramids or similar shapes can be fabricated on a bulk 
substrate. 

In experiments aiming at the investigation of nanoscopic systems we deal frequently 
with such systems deposited on a bulk substrate. However, most often in the theoretical 
description of such experiments the interaction between the nanosystem and the substrate 
is neglected. The main purpose of this work is to find a way to include this and calculate 
its effect on the thermodynamic properties of the Ising model of a ferromagnetic nanoscopic 
system. To achieve this we apply the equilibrium reduced-state (density) operator (ERSO)[^ 
and carry out the calculations within the Gaussian fluctuations approximation (GFA). 

ERSO is the most general equilibrium state operator. It has been derived from the 
generalized Schrodinger variational principle and not on the quantum statistical mechanics, 
i.e. without statistical hypotheses. When the term describing the interaction between the 
system and its environment is neglected ERSO takes the form of the statistical operator of 
the Gibbs canonical distribution. ERSO applied to an exactly solvable microscopic model 
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leads to exact results. 

GFA is a modified version of the high density expansion method as has been proposed in 
refs. 0, S|. GFA is an improvement over the molecular field approximation (MFA) due to 
the self-consistent inclusion of Gaussian fluctuations of this field. The essential new element 
of GFA is the summing up of the partial sums of Feynman diagrams of the same structure 
of recurrent formulae at each stage of the calculations. GFA is based on a classification of 
the Feynman diagrams in terms of where z is the effective number of spins interacting 
with any given spin. Owing to this procedure the theory becomes internally consistent and 
does not lead to unphysical results such as, for example, a complex Curie temperature, {ol 

2. THE MODEL 

In this paper we consider magnetic properties of a nanoscopic pyramid (nanopyramid). 
They are well described by the simple model of localized and ordered spins with the following 
spin-1/2 Hamiltonian 



where / is the coupling parameter and ^ stands for summation over pairs of different simple 
cubic (sc) lattice points. Our considerations are restricted to nearest-neighbor interaction 
only. 

In expression (1) f denotes the two-dimensional position vectors of a spin belonging to a 
given monoatomic layer / = 1,2,3,4 of the pyramid. There are 4 spins in the / = 1 layer, 
16 spins in the I = 2 layer, 36 spins in the / = 3 layer and 64 spins in the / = 4 layer. The 
total number of spins in the pyramid is 120 (see Fig. 1). 

We assume that the nanoscopic pyramid is deposited on a bulk ferromagnetic substrate 
sufficiently well described by the spin-1/2 Ising model for a simple cubic lattice and the 
Hamiltonian 



Here g denotes the two-dimensional position vectors of spin belonging to a given 
monoatomic layer. The summations always run over different sites. The substrate is di- 




(1) 




(2) 
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vided into monoatomic layers parallel to the planes (100) of a sc-lattice. The position of 
each layer is given by the number r = 1, 2, ... 

Let us assume that the interaction of the nanoscopic pyramid with the bulk substrate is 
described by the Heisenberg Hamiltonian of the form 

Hi = ~2^2 Sfi=4Sgr=i. (3) 
fg 

In order to take into account the interaction of the nanoscopic pyramid with the bulk 
substrate we shall apply ERSO^ suitable for description of a physical situation similar to 
the one we are concerned with. In the derivation of this particular form of ERSO we use the 
fact that although the Universe as a whole is in the pure state, its arbitrary multiparticle 
parts are inevitably in mixed states. This is purely a quantum effect following from the 
holistic properties of the quantum theory, formally related to the fact that the Universe 
(according to the quantum cosmology postulates) has one vector of state common for all 
systems and, in the case of interactions among them, we are not able to specify the vector 
of state for individual subsystem (it may not be the case for the interactions of the effective 
field type or classical ones). Such a situation does not occur in the classical description 
because we may know classical trajectories of each particular molecule irrespective of their 
interactions. The holistic features of the quantum theory imply the use of the formalism 
of ERSO in description of multiparticle systems. For this reason the interaction ([3]) of the 
nanoscopic pyramid with the bulk substrate cannot be described by the Ising Hamiltonian. 

In our case ERSO takes the following formQ|: 

d = exp[P{F-H-H'mi (4) 

where F is the free energy and 

H'{(3) = Tts[Hj{1 + K) exp(/5(F, - Hs))] (5) 

is an effective term describing the interaction between the pyramid and the bulk substrate, 
Trs[...] is the partial trace over the substrate states, K is the correlation operator, (3 = 
{ksT)'^, and 

Fs = -^\nTTs[exp{-pHs)] (6) 
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is the free energy of the substrate. Substituting expresssion ([3]) into Eq. ([5]) with the 
assumption that 



I < 5^^^ > I, I < > I < I < > I , /2 < /i , (7) 
and ~ 0, we arrive at 

H\(3) = -\h<SU>Y.^^^i=^ ' (8) 

f 

where 

< SU >= Tr,[^|,=i exp(/?(F, - Hs))] (9) 
is the average bulk substrate spin moment in the layer r = 1. 

3. GAUSSIAN FLUCTUATION OF THE MOLECULAR FIELD 

As a starting point to GFA we choose the following decomposition of the Hamiltonian 

m 

H={H-Hi)+Hi = Ho + Hi, (10) 
where the perturbative part Hi is defined by the tranformation 

H-.Hi = HiS^i^6S^i), (11) 

where 

6S',, = S',, - (12) 
is the fluctuation operator of the z-component of the spin and 

(Sn) = Tr[^f^,e^(^-^-^'(^))], /? = (A^^T)"^ (13) 

and F is the free energy. 

According to the rules of the thermodynamic perturbation expansion we can write 
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(14) 



where 



(...)o = Tr[...rfo] 



(15) 



and 



do = {Tr[exp(-/5(ifo + H'm)]Y' exp(-/3(/7o + H'm)- (16) 

The right-hand side of expression f|T3|) can be expanded into a series with respect to the 
perturbing term Hi. From such an infinite series we now choose a certain partial sum which 
can be represented graphically in the following wayjlo| 



(17) 



where 



il 



+ 



n n n n 

and symbols MM denote the renormalized interaction line, 



f/ f/ 



\ f/ pr n J 



(18) 



(19) 



o o • • . o 



L{Wfi), 



m 



n 



L(iyft) =ln(2cosh(iyn/2)), 



(20) 

(21) 
(22) 



/3 



f'l' 



(23) 
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where we assumed h = 1 and ^fq, denotes summation over nearest neighbours of the U 
spin. 

As a result of calculating the infinite sum (fTSl) and (fT9i) we obtain 



1 

^= / e-'''^^tanh{Wfi + u6Wfi)du, 
V2n J-oo 



(24) 



where 



P / 1 



2 \ 2v^ 



-nV2 



/2 ^(1 - tanh^(iyf,i, + u6Wfi)) + 1261,4(1 - tanh2(K + u6Vr)) 

(25) 

Similarly, using eq. ([9]) for r = 1,2,..., the molecular field K- of the bulk substrate in 
the GFA satisfies the equation 



1/2 



du 



Vr 



\ g gr' / 



(26) 



where 



1 r+00 

(S^) = —= / e-"'/2 tanh(K + m5K)c^m, 
2v27r j-00 



(27) 



= ^(-^ e-"'/2[4(l - tanh2(K + uWr)) 
2 v27r J-oo 



2 

+ (1 - 5r,i)(l - tanh2(\/,_i + ubVr-x)) + (1 - tanh2(K+i + M5K.+i))]rf^)^^^ 

(28) 

and ^PVfi, are the mean Gaussian fluctuations of the molecular fields and l^, re- 
spectively. 

After introducing the following reduced magnitudes. 



equations (24) - fl28|) can be written in a compact form. 



(29) 



-I r+oo 

-= / e-"'/2tanh(iyn + M5Wn)c/M, 

'2tI J-oo 



(30) 
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where 



^'^ = ][Y1 ^f'' + «'^Ml/r=l 1 , (31) 



i V27f 



/ + 00 
e""'/^[l - tanh2(M/n + M'JM/n) + ^(1 - tsmh^Wfn' + uSWfn')) 

(32) 

K = -(4yr + Z/r+1 + (1 - Sr,l)yr-l), (33) 
1 /•+0O 

= ^= / e-"'/^ tanh(K + u5Vr)du, (34) 
V27r 7_oo 



5V; = —{^= / e-"'/^[4(l - tanh2(K + w^K)) + (1 - 5r,i)(l - tanh2(i;_i + -u^K-i)) 

+ (1 - tanh^(K+i + lit^K+OJc^xi)^/^ . 

(35) 

Equations (30)-(35) were solved numerically. In order to carry out the calculations we 
have to specify the magnetic lattices pyramid and substrate as a simple cubic one. Results 
are presented graphically. 

Fig. 2 shows the spatial distributions of magnetization X4 in the I — 4 layer of the 
nanoscopic pyramid for temperature t = 4.0 and a = 0.5, b = 1.0 in both MFA and GFA. 

Figs. 3-5 present spatial distributions of Gaussian fluctuations of molecular field 6Wfi for 
/ = 2, 3, 4 respectively, temperature t = 4.0 and a = 0.5, b = 1.0. In the tip layer adjacent to 
the substrate the magnitude of fluctuations is largest in the second row of atoms, counting 
from the outside. Comparing Fig. 4 and 5 we see that maximum fluctuations occur along 
the edges of the pyramid. Smaller fluctuations of magnetic moments in the center of the 
base of the apex is an effect of the substrate. 

Fig. 6 shows the mean value of magnetization 

Xi^^^ , (36) 
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TABLE I: Curie temperature tc (in relative units) for the nanoscopic pyramid. 



Qj 


b 


MFA 




GFA 





cxL uiuicii y 








0.1 


0.5 


4.93 




1.88 


0.1 


1.0 


5.99 




4.13 


0.1 


2.0 


11.95 




9.30 


0.5 


1.0 


5.99 




4.13 


0.5 


2.0 


11.95 




9.27 



in each monoatomic layer of pyramid Z = 1, 2, 3, 4 in GFA as a function of temperature t for 
a = 0, 6 = 1 and a = 0.5, b = 1.0. 

Fig. 7 presents temperature dependence of the mean pyramid magnetization 



X = (4Xi + I6X2 + 36X3 + 64X4)7120 (37) 

obtained in GFA for a = and a = 0.5 and b = 0.5, 1.0, 2.0. Fig. 8 shows magnetization of 
individual layers of the pyramid and the substrate. The influence of the pyramid extends 
several layers into the substrate. Similar conclusion was reached in an earlier study 121] of 
the interaction between a magnetic nanotip and a magnetic surface using the tight-binding 
model. The tight-binding calculation for a Fe tip on a Fe surface shows that only the first 
four layers of the tip support are affected. The inclusion of fluctuations in our work increase 
the depth of this influence as t — > tc, see Fig. 8. 

Fig. 9 presents the mean Gaussian fluctuation of molecular field 

SW, = (38) 

and 6Vr for pyramid and the substrate respectively. 

Finally, Fig. 10 shows Gaussian fluctuation of the molecular field as a function of tem- 
perature. 

In Table 1 the Curie temperature tc of the ferromagnetic nanoscopic pyramid is given in 
relative units for different values of parameters a and b obtained in MFA and GFA. 
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We can see that for a nanopyramid the ratio of the Curie temperature in GFA to that 
obtained in MFA is 0.783. For an infinite simple cubic lattice this ratio is equal to 0.856 

n 

and for a square lattice Ising monolayer it has a value of 0.799. [8| The ratio of the Curie 
temperature of a monolayer obtained in GFA to the exact result is 1.409. The same quantity 
for the sc bulk Ising ferromagnet is approximately 1.138. 

4. CONCLUSIONS 

The results we obtained lead us to the following conclusions: 

• The interation of a ferromagnetic nanoscopic pyramid with its bulk ferromagnetic 
substrate may have essential influence on the properties of the ferromagnetic pyramid, 

• The distribution of Gaussian fluctuations is highly nonuniform, 

• limj^o SWfi = 0, as expected, 

• The maximum of molecular field fluctuations SWfi is reached at the Curie temperature 
tc, as expected, 

• The Curie temperature tc of the pyramid strongly depends on the Curie temperature 
of the bulk substrate, 

• We believe that ERSO (jl]) may be successfully applied in studies of the influence of the 
bulk substrate on the thermodynamic properties of nonmagnetic systems and other 
more complex nanoscopic systems. 

m 



The development of the magnetic exchange force microscopy (MExFM) 13|] may enable 
studies of the spatial distribution of magnetization in nanopyramides deposited on nonmag- 
netic surface. Such magnetic particles are of interest in magnetic recording techniques. 

It would be useful to study the temperature dependence of exchange force between the 
nanotip and a well-characterized magnetic surface of a material with Curie temperature 
higher than of the tip. On cooling the system one should observe the magnetic transition 
of the tip. The force of the interaction between the tip and the surface is measured as the 
shift of frequency of oscillations of the cantilever. As T decreases, magnetization of the tip 
increases and so does the exchange interaction with the surface, resulting in a frequency 
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shift. If the ratio of exchange energies /1//2 is large, i.e. when exchange interaction between 
atoms within the tip is much larger than exchange between the bottom layer of the tip 
and the first layer of the substrate, or when the substrate is nonmagnetic, fluctuations may 
lead to a long tail of magnetization as a function of temperature, see Fig. 7. In this case 
magnetization may saturate at temperatures much smaller than T^. 

The measurement of the frequency shift as a function of temperature at fixed position 
above the surface may provide a good estimate of possible magnetic fluctuation effects. The 
crossover temperature between the tail, see 6 = 2 curve in Fig. 7, and the saturated part at 
low T detected via the oscillation frequency shift should not depend on position above the 
surface. Temperature scans at different surface positions would result in frequency shifts of 
different magnitudes occurring at approximately the same temperature. 

Recent attempts to map single spins on the surface of antiferromagnetic insulator NiO 
were only partially successful. The experiment of Kaiser et al. 13| with Fe tip performed at 
about lOK in a magnetic field of 5 T showed clear picture of antiferromagnetic structure of 
the (001) surface of nickel oxide. It cannot be ruled, however, that the applied magnetic field 
induces structural changes in NiO. The measurements of Schmid et al. 14\ conducted at room 
temperature, using Co and NiO tips did not did not show any spin contrast. The ordering 
temperature of NiO is 525 K. The authors point out that the magnetic ordering temperature 
of the tip may be significantly reduced relative to the bulk value. Our calculation shows that 
it is also possible that the spin ordering of the tip is present jDut is very weak. Performing 
this experiment at lower temperatures, as authors of Ref. [xJ] suggest, may yield better 
results. 

Atomically sharp magnetic tips are crucial in achieving lateral high spatial resolution. 
Therefore the knowledge of spatial distribution of fluctuations and its temperature depen- 
dence might improve the calibration of the STM device and help in estimating experimental 
error of the MExFM method. For small distances between the tip and the surface it may be 
necessary to include interaction of more than one atom of the tip or of the surface. 151] In this 
case our model with nonuniform magnetization of the tip may provide useful insight. The 
temperature dependence of mean value of magnetization of the individual tip layers shown 
in Fig. 6 shows that the magnetization in the second layer can be significantly larger that 
the magnetization of the apex. This implies an increased range of tip-surface separation 
where the interaction with atoms in the next layer should be taken into account. 
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The ability to control magnetic properties of nanoclusters[16!] is one of the central prob- 
lems of nanotechnology. Arrays of nanoclusters are also intensively studied, fl?!. [isl The 
theory presented here may be applied to various geometrical shapes used in experiments, 2I 
e.g. chains of particles, striped and cylindrical nanowires, nanodots, nanojunctions, surface 
steps. We intend to carry out calculations in a more realistic model of the substrate surface 
and the interface between the pyramid and the substrate. 

Finally let us note that ERSO (jl]) couples the nanoscopic pyramid with the substrate. 
This means that the system does not have a finite number of spins. Therefore the use of 
e.g. Monte Carlo simulations would have little justification. 
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FIG. 1: Schematic diagram of ferromagnetic nanoscopic pyramid deposited on a ferromagnetic 
bulk substrate. 
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FIG. 3: Spatial distribution of Gaussian fluctuations of the molecular field in pyramid layer 1 = 2 
at temperature t = 4.0 and a = 0.5, 6=1. 
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FIG. 4: Spatial distribution of Gaussian fluctuations of the molecular field in pyramid layer I = 3 
at temperature t = 4.0 and a = 0.5, 6=1. 
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FIG. 5: Spatial distribution of Gaussian fluctuations of the molecular field in pyramid layer 1 = 4 
at temperature t = 4.0 and a = 0.5, 6=1. 
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FIG. 6: Mean value of magnetization, Xi = ^f^^fi/Xlf ™ each monoatomic layer of pyramid, 
I = 1,2, 3, 4, in GFA as a function of temperature i for a = 0.5, b = 1.0. 
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FIG. 7: Temperature dependence of the mean pyramid magnetization X = (4Xi + I6X2 + 86X3 + 
64X4)7120 for a = 0.5 and b = 0.5, 1,2. 
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FIG. 8: Magnetization in individual layers of the pyramid Xi and the substrate yr for a = 0.5, 
b= 1.0 and t = 4.0,4.12 
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FIG. 9: Gaussian fluctuation of the molecular field in individual layers of the pyramid and the 
substrate for a = 0.5, b = 1.0 and t = 4.0, 4.12. 
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